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We compute the exact single-particle time-resolved dynamics of electronic Mach-Zehnder inter¬ 
ferometers based on Landau edge-states transport, and assess the effect of the spatial localization 
of carriers on the interference pattern. The exact carrier dynamics is obtained by solving numeri¬ 
cally the time-dependent Schrodinger equation with a suitable 2D potential profile reproducing the 
interferometer design. An external magnetic field, driving the system to the quantum Hall regime 
with filling factor one, is included. The injected carriers are represented by a superposition of edge 
states and their interference pattern reproduces the results of [Y. Ji et al. , Nature 422, 415 (2003)]. 
By tuning the system towards different regimes, we find two additional features in the transmission 
spectra, both related to carrier localization, namely a damping of the Aharonov-Bohm oscillations 
with increasing difference in the arms length, and an increased mean transmission that we trace to 
the energy-dependent transmittance of quantum point contacts. Finally, we present an analytical 
model, also accounting for the finite spatial dispersion of the carriers, able to reproduce the above 
effects. 
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I. INTRODUCTION 

Edge states (ESs) are chiral one-dimensional conduc¬ 
tive channels which arise in a 2D electron gas (2DEG) 
in the Integer Quantum Hall regime (IQHE) 1,2 . They 
are essentially immune to backscattering and are char¬ 
acterized by very long coherence lengths 3 . Besides 
their remarkable interest for basic solid-state physics 
and coherent electronic devices, they are ideal candi¬ 
dates for both the so-called “electron quantum optics” 4 
and “flying qubit” implementations of quantum com¬ 
puting architectures 5,6 . 

Indeed, several attempts to demonstrate single-qubit 
gate operations and electronic interference in coupled 
quantum wires were hampered by scattering and de¬ 
coherence processes 7 9 . On the other hand, experi¬ 
mental realizations of electronic interferometry based 
on edge-channel transport seem more mature, to the 
point of demonstrating not only single-electron 10,11 
but also two-electron interference 12 . Also, “which- 
path” detectors or quantum erasers 13-17 have been 
implemented and the formation of quantum entan¬ 
glement between indistinguishable particles has been 
demonstrated 12,18-21 . 

While specific features of electronic Mach-Zehnder In¬ 
terferometers (MZIs) 22 are still under investigation 
and involve possibly many-particle effects 21,23,24 , their 
basic functioning can be explained in terms of station¬ 
ary electron waves interference, like their more com¬ 
mon optical counterpart. However, contrary to pho¬ 
tons, the time that a conduction electron takes to cross 
the interferometer can be comparable to switching 
times of typical microelectronic devices, and its spatial 
localization can be much less than the dimensions of a 
single arm of the interferometer. Thus, understanding 
the detailed quantum dynamics of the carriers, beyond 
the simple plane wave model is critical for the design 
of novel devices based on ES transport. This is even 
more relevant if the device is intended to process the 


information encoded in a single particle at a time, a 
task which requires a high resolution both in time and 
space. 

In fact, the injection of single-electron excita¬ 
tions that propagate along ESs has been recently 
demonstrated 11,25,26 : with a proper time modulation 
of the injecting pulse, they consists of Lorentzian- 
shaped excitations above the Fermi sea, termed levi- 
tons 27 29 . However, leviton-based MZI has not been 
realized experimentally so far. The majority of the 
literature on electron MZIs in ESs deals with elec¬ 
tronic currents 10,12,13,17,22 , and the prevalent theoret¬ 
ical models for the transport are based upon delo¬ 
calized scattering states 14-16,18,30,31 , with few notable 
exceptions 32,33 . 

In this work we address, both from numerical and an¬ 
alytical points of view, the interference properties of 
a Mach-Zehnder device based on ES channels tuned 
to filling factor v=l and quantum point contacts used 
as splitting elements, when the electron travelling in¬ 
side it is strongly localized. We first use a numeri¬ 
cal approach, which, unlike other recent works 34,35 , is 
based upon the direct solution of the effective-mass 
time-dependent Schrodinger equation, in presence of 
an external magnetic field. The electrons travelling 
inside the device are localized wave packets (WPs) of 
ESs, and are propagated by a generalized split-step 
Fourier method 36 38 . The effects of the WP size on 
the transport process is analyzed. Then, we develop 
an analytical model also accounting for the finite spa¬ 
tial dispersion of the carriers. The transmission coef¬ 
ficient of the device subject to Aharonov-Bohm (AB) 
oscillations 39 is obtained as a function of the magnetic 
field and of the geometrical parameters of the inter¬ 
fering paths, and the results are compared with the 
numerical simulations. Specifically, in Section II, we 
describe the initial electronic wave function and the 
physical device used in our simulations. In Section 
III we give some details on the numerical algorithm 
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that we adopt and we describe the results of the nu¬ 
merical simulations. Section IV is devoted to the de¬ 
velopment of an analytical model for the transport, 
which takes into account the energy dispersion of the 
wave packet. A comparison between the predictions 
of the latter model and the numerical simulations is 
then presented in Section V. Finally, the conclusions 
are drawn in Section VI. 


II. PHYSICAL SYSTEM 

We consider a conduction-band electron in a 2DEG 
on the x-y plane, with charge — e and effective mass 
m*. A uniform magnetic field B = (0,0,5) is ap¬ 
plied along the z direction and a non-uniform electric 
field induces a local potential energy V(x,y). The 
latter will represent the field generated by a polarized 
metallic gate pattern above the 2DEG that will define 
energetically forbidden regions. The generic Hamilto¬ 
nian H = {—ihV + eA) 2 /(2m*) + V can be rewritten 
in a more explicit form by substituting the magnetic 
vector potential A with its expression in the Landau 
gauge A = (0, 5x, 0). The 2D Hamiltonian for the 
electron in the 2DEG then reads 

__ h 2 d 2 . h eBx d 

2m* dx 2 2m* dy 2 m* dy 

+ v { X ,y). ( 1 ) 

We adopt a single-particle approach, thus neglecting 
the interaction with other free electrons of the de¬ 
vice. Our time-dependent numerical simulations and 
analytical model are based on Eq. (1), as we explain 
in the following. Before that, we must recall briefly 
the derivation of Landau states and of corresponding 
ESs 40 . 

Let us consider a region where the electric poten¬ 
tial is invariant along y, i.e. V(x,y) = V(x). In 
the Landau gauge, the Hamiltonian (1) shows an ex¬ 
plicit translational symmetry along the y direction, 
and the quantum evolution of the particle can be fac¬ 
torized along the two axes. In this case, the eigen¬ 
states of H have the form T(x,t/) = Lp{x)e lky . In 
fact, the ^-dependent part of the wave function (WF) 
is a plane wave, while the x-dependent part (which 
depends also on the wavevector fc) is a solution of 
P n ,k( x ) = F n (k)ip ni k(x), with the following ID 
effective Hamiltonian 

+\™?u 2 c (x~xq) 2 + v (x), ( 2 ) 

where 

hk 

= ( 3 ) 

and where ui c = is the cyclotron frequency. Note 
that the parameter xo representing the center of the 
effective parabolic confinement along x depends on fc, 
i.e. the wavevector in the y direction. The discrete 
energy levels E n (k ) associated with are the so- 

called Landau levels (LLs). 


If V(x) = 0, the system eigenfunctions are the well- 
known localised Landau states. However, if V(x) is 
a step-like function, defining a sub-region in which 
the electrons are confined, the Landau states with xo 
close to the edge have higher energy, and show a fi¬ 
nite dispersion in k. They become edge states , which 
are delocalized WFs associated with a net probabil¬ 
ity density flux, and which can act as ID conductive 
channels. 

Since we want to model a carrier as a propagating 
wave packet, beyond the delocalized scattering state 
model, we need to construct such WF as a proper 
combination of ESs, rather than considering a single 
ES. This is similar to representing a free electron by a 
minimum uncertainty wave packet rather than a plane 
wave, thus including a finite uncertainty in both its 
position and its kinetic energy. However, in our case 
the magnetic field couples the two directions and some 
care has to be taken, since the position along x and 
the velocity along y are not independent. 

We suppose that the injected electron lies in the first 
LL and we represent it as a linear combination of ESs 
with n — 1, with Gaussian weights. We choose a 
Gaussian shape for our initial WP since it represents 
a quantum particle with the minimum quantum un¬ 
certainty in both its position and momentum, and be¬ 
cause its longitudinal spreading with the time is lim¬ 
ited. This allows us to follow clearly the dynamics 
of the localized WP without substantial high-energy 
components forerunning the center of the WP and 
leading to numerical instabilities. This choice also 
allows the derivation of the approximated analytical 
model presented in Sect. IV. However, the qualita¬ 
tive results of our simulations do not depend on the 
specific shape of the initial WP. In fact, although we 
address the experimental regime of Ref. 10 and not 
the propagation of levitons, as in Ref. 25, we tested 
different shapes for the initial electron wave function: 
due to the superposition principle the final state will 
be the composition of the different scattering states 
taken with their initial weight. 41 Specifically, the wave 
function at the initial time will be 

*o(x,v) = j dkF(k)<p hk (x)e iky , (4) 

where F(k) is the Fourier transform of a ID Gaussian 
function along the y axis: 


F{k) = LJ d y e~ iky 


( y-vo ) 2 , 

e~ 4 CT 2 e + lkoy 


\Z~2l TO 2 


( 5 ) 


27T 3 


e -a 2 (k-k 0 ) 2 e ~ikyo 


This leads to a localized WF along both directions: 
in t/, the Gaussian envelope gives a finite extension 
around the central position 7 / 0 , in x, the functions 
(pi,k{x) are always localized around xo(fc), and the 
wave vector k is, in turn, localized around fco by our 
choice of F(k). With the initial condition Eq. (4), our 
simulations will be able to take into account the effects 
of the size a of the wave packet (WP) on the inter¬ 
ference phenomena. In order to use Eq. (4), we will 
be careful to initialize the WP 4 /q where the potential 
V (x, 7 /) does not depend on y. 
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Figure 1. (Color online) Device geometry and square mod¬ 
ulus of the WF at t — 0 (initial state, with a = 40 nm), 
t = 2.5 ps (at the first QPC), t — 7 ps (along the two 
arms), t = 12.5 ps (at the second QPC). In this example 
B — 5 T and — —65 nm. The resulting transmission is 
T = 0.599. 


The specific form of the ESs depends on the the shape 
of the potential barrier V(x), and in general there 
is no closed-form expression for them. In order to 
have a realistic model of the smooth edges between the 
allowed 2DEG region and the depleted one, defined 
e.g. with the split-gate technique, we take 

V(x) = VoTr(x), (6) 

where J r T (x) = (exp(rx)+ 1) _1 is a Fermi distribution 
with a “broadening” parameter r, and we compute 
numerically the corresponding <p n ,k(x ) states. 

The shape of the depleted regions of our 2DEG is 
chosen in order to mimic, with the ES channels, the 
Mach-Zehnder interferometer (MZI) of Ref. 10. This 
is depicted in Fig. 1, where the dark and light regions 
represent V(x,y) = 0 and V(x,y) = Vo, respectively. 
Two narrowings, where two areas with high poten¬ 
tial are close to each other, form two quantum point 
contacts (QPCs), both with a square area 42 of S 2 : 
their dimension S is tuned in order to give a trans¬ 
missivity of 50% in the two output channels. The 
localized electron is injected from the top of Fig. 1 
(where the initial WF is centered around x = — 20 nm 
and y = 700 nm), in the first LL, whose ESs follow 
the boundary between the high- and low-potential re¬ 
gions. After the first QPC, each ES is split into two 
parts that constitute the two arms (“left” and “right”) 
of the interferometer. They are rejoined at the sec¬ 
ond QPC. As a consequence, the two components of 
the WP that follow the two arms, interfere at the sec¬ 
ond QPC. Indeed, two outputs are available there to 
the electron: one part is reflected inside the MZI, and 
then absorbed by a contact (modeled as an absorbing 
imaginary potential V a bs, as detailed later), the other 


part is transmitted towards the bottom of the device. 
The norm of the latter part gives a measure of the 
transmission coefficient T of the device. 

The geometrical parameters of the MZI are indicated 
in Fig. 1. Specifically, the “right” arm and the “left” 
arm have two horizontal sections each, of length Lq 
and Lq T (2.l, respectively. The vertical sections of 
the two arms have obviously the same length, Hq. 
The absorbing lead is modeled by a purely imaginary 
potential 38 ’ 43 


Vabs( x iV) ~ abs 


0 r( X Xc)T rix X c S c ) 


cosh 2 ((y - y c ) 2 /d 2 ) 


( 7 ) 


where x c and y c are the coordinates of the contact, 
while S c and d are, respectively, the contact extension 
along the x and y axes. By changing the parameters 
of the device, such as (which controls the relative 
length of the two arms and the area) or the magnetic 
field H, we observe AB oscillations in the transmission 
coefficient T. With the introduction of the absorbing 
lead V a bs , the calculation of T is straightforward, since 
it is the norm of the final wave function after the ab¬ 
sorption process. Differently to other works 44 focusing 
on the dynamics of carriers moving along Hall ESs, our 
confining potential V does not depend on time. 

The “center” x 0 (&o) of the initial WP (t = 0) is fixed at 
a distance A from the inflection point of the potential 
barrier V(x), and consequently fco can be calculated 
by Eq. (3). The resulting local bandstructure of the 
1st LL around k = fco can be determined numerically, 
and from a parabolic fit we get the values of the pa¬ 
rameters Eq, m B and fci in its expression: 


E^k) 


h 2 

2 m* B 


(k-hf + Eo. 


(8) 


The parameters fco and fci are gauge-dependent quan¬ 
tities, since they depend by the origin of the coordi¬ 
nate system: however, with a proper choice of this 
degree of freedom and of the energy zero, we can set 
fci = 0 and Eq = 0 without any physical change in the 
description of the system, in order to get a simpler de¬ 
scription of the system. Now fco is directly related to 
the group velocity v g = hko/m* B of the WP, which 
behaves like a free particle of mass m B in ID. 

In the simulations we use for the electron the effective 
mass of GaAs m* = 0.067 m e . The parameters 
of the potential and of the absorbing contact are 
Vo = 10 eV, t -1 = 3 nm, S c = 100 nm, d = 30 nm, 
Va bs = — 08.7 eV; the constructive parameters of 
the device are S = 32.2 nm, Lq = 350 nm and 
H 0 = 600 nm; the parameters of the WP are 
A = 20 nm and a = 20,40,60 nm 45 , while the 
magnetic field is kept around the value of B = 5 T. 


III. NUMERICAL SIMULATIONS 

The time evolution of the WP is realized with the 
split-step Fourier method 38,43,46 . In summary, the 
evolution operator U(St) = e~i 5t ' H is applied N times 
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to the initial wave function 4/(x, y, 0), each leading to 
the evolution of a short time step St: 

t)\ t =NSt = [U{5t)] N y(x,y, 0). (9) 


T2(Px) = 


h 2 d 2 

2m* dx 2 


(ii) 


The Hamiltonian of Eq. (1) can be written as H = 
Ti(x,p y ) +f 2 {p x ) + V(x,y), with 


^1 Py ) — 


h 2 d 2 
2m* dy 2 


cBx d e 2 B 2 9 

ih 1 


m* dy 2m 


<(10) 


By applying the Trotter-Suzuky factorization and the 
split-step Fourier method, the operator U(St) can be 
approximated with 




[U(St)r = 


,-kst-(y+v abs )^ i e -kst-t 1 &0--i e -i5t-f 2 & 


N 


e * dt 2 


( 12 ) 


where and ^x(i) denote, respectively, the di¬ 

rect and inverse Fourier transform with respect to the 
variable x(y). By using the above expression, the nu¬ 
merical solution of Eq. (9) is reduced to an alternating 
application of discrete Fourier transforms and array 
multiplications, since each operator acts only in its 
diagonal representation. 

Within this scheme, we can reproduce the exact dy¬ 
namics of the electronic WF inside the device. As 
an example, the evolution of the WP at four different 
times is shown in Fig. 1 for the specific case described 
in the caption. The simulations have been performed 
with a time-step St = 0.1 fs. 

As expected 10,12,17 , we obtain AB oscillations of the 
transmission T (see Figure 2) with respect to both 
the variation of the area (which is controlled through 
cll) and the variation of the magnetic field B. This is 
consistent with the results of Refs. 10 and 12, where 
the tailoring of the length of one of the two paths is 
achieved through a modulation gate. 

However, we observe two extra features, namely (a) 
a damping of the AB oscillations with ac, and (b) an 
increase of the average T as the initial spatial local- 
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Figure 2. (Color online) Aharonov-Bohm oscillations in 
the transmission coefficient T as a function of the asym¬ 
metry parameter of the left arm aL and of the magnetic 
field B (numerical simulations for a wave packet with 
a = 40 nm). 


ization a of the WP is increased, with a consequent 
decrease of the visibility 10,12,13,22 


VMZI = 


Imax I min 
Imax T I min 



(13) 


of the AB oscillations (of the transmitted current /), 
which is always smaller with respect to the ideal case 
vmzi = 1. The two issues above are addressed in the 
following. 

(a) The AB oscillations as a function of are mod¬ 
ulated by a Gaussian envelope, whose extension is di¬ 
rectly correlated with the size a of the initial WP. A 
fitting of T vs a is reported in Fig. 3, for three values 
of a (see caption). This effect has a simple physical 
explanation, already advanced by Haack et al. in 32 
for a Lorentzian WP. Indeed, when the asymmetry be¬ 
tween the two arms is large with respect to the size of 
the WP, the two parts of the WF arrive at the second 
QPC at different times, and do not interfere. In this 
case, each part is transmitted with a probability of 
50%, ending up with a total transmission of T ^ 0.5. 
In general, the larger the time offset at which the cen¬ 
ters of the two WPs reach the second QPC, the less 
effective is the quantum interference. Therefore, we 
expect a saturation value of ~ 0.5 for T(cll ), and an 
oscillation amplitude of ^ 0.5. However, our numeri¬ 
cal simulations show a different trend, i.e. issue (b). 

(b) For smaller values of cr, the average (or saturation 
value with cll ) of the transmission T is higher and the 
amplitude of the oscillations is lower (Fig. 3). These 
effects are due to the energy-dependent features of the 
scattering process at the QPCs, which are included 
automatically in the numerical simulations based on 
the direct solution of the Schrodinger equation. In 
fact, contrary to a delocalized scattering-state model, 
where the particle is represented by a single-energy 
state, our WP is composed by different ESs, as given 
in Eq. (4), each with a slightly different energy. As 
a consequence, the two parts of the WP, transmit¬ 
ted and reflected by the QPC, have different spectral 
weights, which depend on a. Only in the limit of 
<r oo, the WP is split into two identical parts that 
give an ideal 50% splitting with T = 0.5. This effect 
should be detectable in noise spectra of two-particle 
scattering, as proposed in Ref. 47 . In order to have 
a better physical insight and to give a quantitative 
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Figure 3. (Color online) Transmission coefficient T as a 
function of the asymmetry parameter of the left arm aL 
(top curve: a = 20 nm, middle curve: a = 40 nm, bottom 
curve: a = 60 nm). The curves are the fits of the numerical 
simulation data. It turns out that the best fit is obtained 
with the function T(cll) — c\ + C 2 • exp[— c^{aL — C 4 ) 2 ] • 
cos(c 5 aL + C6), where the ci are the fit parameters. Inset: 
Detail of the curve cr — 40 nm with numerical simulation 
data-points. All the simulations are performed at B = 5 T. 


assessment of this effect, we will include the energy- 
dependent transmissivity of the QPCs in the analyti¬ 
cal model presented in the following section. 


IV. THEORETICAL MODEL FOR THE 
INTERFERENCE 

Following Weisz et al 17 , we will represent a single ES 
as a ID plane wave, and we will describe the scattering 
process and the transport inside the MZI with the 
scattering matrix formalism 2 . Of course, the effects of 
the scattering on the WP can be determined from the 
linear superposition of the scattering of single plane 
waves. 

The device can be divided into three regions - I, II 
and III - as indicated in Fig. 4, where the plane-wave 
states are given by: 

1 . I *M*0 >-(J) (14) 

II. | L(k)) = ( J ) \R(k)) = ( J ) , (15) 

III. |A2(fc)) = ( J ) KM*)) = ( J ) • (16) 

The initial wave packet is prepared in the state 

\*i) = J dkF(k)\ Dl (k)), (17) 

and the scatterings at the two QPCs are described by 
the scattering matrices 

*=(**)• (i8) 


I A) I 


II 


l^> 

- w 


! f i \ u *) 

' 1 - 


_1 - L._ 

III 

|A> | 


Figure 4. Scattering states inside the device. Regions I, II 
and III are separated with dashed lines. 


with j = 1 (j = 2) labeling the first (second) QPC, 
and where rj and tj are the transmission and reflection 
amplitudes. After the first QPC, the WF is given by: 

I*//) = S 1 1*/) 

= I dk (. F L (k ) I L{k)) + F R {k) \R(k))) , (19) 

where F^k) = r\{k)F(k) and Fji(k) = iti(k)F(k). 
While travelling in the two arms of the MZI inside 
region II, the WF acquires a different phase in each 
arm, which is given by 2,48,49 : 

<Ph = J{p — qA) ■ ds = £h+ <t>h, ( 20 ) 

where h = L {h = R) labels the left (right) arm and 
the integral is a line integral along the arm h. The 
first contribution, the dynamic phase Q*, is due to the 
canonical momentum p, while the second contribu¬ 
tion, the magnetic phase accounts for the vector 
potential A. The matrix describing the phase acquired 
is given by 


P = 


e^ L 0 
0 e^ R 


( 21 ) 


and therefore the final WF in region III is given by 

\*in) = S 2 P\Vn) 

= IdkF(k) (t\D 2 (k))+r\U 2 (k))), (22) 


where t = — txt^e 1 ^ 11 and r = ir^e 1 ^ 1 ^ — 

it\r 2 e lipR . It is easy to see that the transmitted part 
of the WP is given by 


|*m 


>-(sP 


- dk'\D 2 (k'))(D 2 (k')\ |* m >, (23) 


where (D 2 (k')\D 2 (k)} = 2irS(k — k') : and therefore 
T-(9j n \*ni), 

= 27 r J dk \F(k)\ 2 \t(k)\ 2 . 


(24) 
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If the scattering amplitudes rj and tj are energy- 
independent, Eq. (24) can be integrated analytically, 
to give: 

~ A l 2 

T = T 0 — Tie 8^2 cos (A(/> — fc 0 A/), (25) 

where T 0 = |rir 2 | 2 + \tit 2 \ 2 and T x = 2r 1 r 2 t 1 t 2 . T 0 
is the mean value of the oscillations, while T\ is their 
maximum amplitude. In the previous formula, A l ~ 
2a l is the length difference between the two arms of 
the interferometer, while A0 ~ j.BHo(aL 0 + 2 L) = 
7r^- is the well-known Aharonov-Bohm phase 2,40,48 
(&b is the flux of the magnetic field through the area 
of the MZI and 4>o = ^ is the magnetic flux quan¬ 
tum) 50 . 

V. DISCUSSION 

As one can see, Eq. (25) provides the correct trend for 
the results of the simulations, with a sinusoidal oscil¬ 
lation enveloped by a Gaussian and with a fixed offset 
in the transmission. Indeed, by substituting the ex¬ 
pressions for Aand A l into Eq. (25), we see that the 
transmission coefficient should oscillate as a function 
of aL like cos(k e aL + tfo), where k e = %BH — 2ko ~ 
4.80 • 10 9 m -1 , which is consistent within 2% with the 
values of k e obtained by fitting the simulation data 
(see Table I). For what concerns the amplitude and 
the saturation value of T, if we suppose that the QPCs 
have a perfect half-reflecting behavior for all incoming 
energies - that is tj = r j = ^ for k - then our 
model predicts a saturation value To = 0.5 and an os¬ 
cillation amplitude T\ = 0.5, as we previously stated. 
This prediction, however, does not match with the 
numerical simulations, where we observe higher values 
for T 0 and lower values for T\ (see Table I) . The model 
also predicts that the oscillations of T{cll) should be 
modulated by a Gaussian with a standard deviation 
E = cr, i.e. with the same spatial dispersion of the 
initial wave packet. However, the actual values of E 
obtained from numerical fitting the data of Fig. 3 are 
bigger than a by more than a factor of two in the case 
with smaller cr. As cr increases, the values of E (always 
obtained by fitting the results of the simulations) get 
closer to cr (see Table I). Indeed, a better agreement 
can be obtained by enhancing our model. 

This simplified model is valid under certain approx¬ 
imations, that are: (i) the scattering process can be 
treated as a quasi-ID problem; (ii) all the edge chan¬ 
nels involved in the transmission follow the same path, 
i.e. we can consider the area enclosed by the interfer¬ 
ing paths and the difference in length A l as they were 
independent on fc; (iii) the reflection and transmis¬ 
sion amplitudes rj and tj of the QPCs are energy- 
independent (that means also /c-independent). While 
assumptions (i) and (ii) are in general verified, the ap¬ 
proximation (iii) must be dropped in order to repro¬ 
duce the results of the simulations. Indeed, our re¬ 
sults can be explained considering that the scattering 
amplitudes of the QPCs are not energy-independent, 
as it has been already pointed out in the literature 19 . 
From a physical point of view, the edge states of higher 


energy are closer to the barrier, and therefore they 
should tunnel easier through the QPCs. This predic¬ 
tion can be easily verified, since the energy-dependent 
scattering process should produce different values in 
the weights Fi,{k) and Fn(k) of the WP after the first 
QPC. This can be linked directly to the differences 
in the values of To and T\ with respect to the ideal 
case, as we will see briefly (a model for the correct 
estimation of E, Tq and T\ is given in Appendix A). 


Initial WP 

T(a L ) (fit) 

a (nm) 

To 

Ti 

E (nm) 

k e (nm 1 ) 

20.0 

0.7371 

0.2403 

40.98 

4.86 

40.0 

0.6330 

0.3650 

54.54 

4.86 

60.0 

0.5778 

0.4215 

70.99 

4.86 


Table I. Data of the initial WP (a is the standard devi¬ 
ation of Eq. (5)) and the corresponding parameters for 
the transmission coefficient T(cll) profile. Data were fit¬ 
ted with the expression T(ol) = To — Ti exp(— \(gll — 
c) 2 /E 2 ) cos(/c e OL + tpo), taken from Eq. (25). 


The amplitudes t(k) and r(k) (which are equal for 
both QPCs, due to their identical shape) can be calcu¬ 
lated numerically from |T//). In Eq. (19), the squared 
modulus of the weights \Fi,(k)\ 2 and \Fn(k)\ 2 has been 
calculated from a simulation on the wave packet with 
a = 40 nm by projecting |T/j) over the local eigen¬ 
states | L(k)) and |i?(fc)) 51 . Then, using a phenomeno¬ 
logical approach, the corresponding data have been 
fitted with the functions /^(fc) = F ar {k — k r ) \F(k)\ 2 
and /#(&) = F at {k — k t ) |T(/c)| 2 , thus giving a Fermi- 
distribution expression 52 for \r(k)\ 2 and \t(k)\ 2 . The 
results of the fits shown in Fig. 5 are in good agree¬ 
ment with the numerical simulations and are also con¬ 
sistent with the constraint \r(k)\ 2 + \t(k)\ 2 = 1. From 
the same figure, we also deduce that the QPCs are 
more transparent at higher energies (i.e., lower values 
of fc), as we predicted. 

With the hypothesis that t(k) and r(k) are real 17 , we 
can insert the transmission and reflection amplitudes 
calculated numerically into Eq. (24) and find the ex¬ 
act transmission coefficient T. The results for a WP 
of cr = 20 nm are shown in Fig. 6: here, the simplified 
transmission coefficient of Eq. (25) (where t and r are 
energy-independent and both equal to 1 /a/ 2) is com¬ 
pared with the exact energy-dependent calculation of 
Eq. (24), which reproduces the result of numerical 
simulations (see Fig. 3 for comparison). Then, as we 
see, the energy dependence of the scattering ampli¬ 
tudes t(k) and r(k) is able to account for the values 
of To and T\ and also for the broadening of E (the re¬ 
gion in which we observe interference) with respect to 
the original spatial extension cr of the WP. Therefore, 
taking into account these corrections, it is still possi¬ 
ble to use the simplified model of Eq. (25) to describe 
the interference, provided that one uses the energy- 
dependent values for Tq, T\ and E, that we name Tq, 
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Figure 5. (Color online) Weights |Tr(/c)| 2 and |Tl(/c)| 2 
for the transmitted and reflected part of |T//) obtained 
from numerical simulations (<r = 40 nm); the continuous 
lines are the fits with the functions fii(k) and /^(fc), the 
dashed line is \F(k)\ 2 (which is peaked on k = ko). In¬ 
set: Squared modulus of transmission and reflection am¬ 
plitudes t(k) and r(k) for a single QPC calculated from 
numerical fits. 


Ti and E. They can be calculated from our model as: 

T 0 = J dk \F(k)\ 2 (|nr 2 | 2 + |tit 2 | 2 ) , (26) 

Ti — 27r J dk\F(k)\ 2 2rir2tit2, (27) 

(see Appendix A for E and for the details of the deriva¬ 
tion of To and Ti). In summary, we obtained an an¬ 
alytical formula for the visibility vmzi of the device, 
which is 

vmzi = A (28) 

-*0 

Indeed, we see that the device has a reduced visibil¬ 
ity with respect to the expected ideal value vmzi = 
Ti/Tq = 1 , due to phase-averaging effects between dif¬ 
ferent energies involved in the WP (and not because of 
decoherence effects), as some experimental works al¬ 
ready pointed out 10 . Interestingly, these remarkable 
effects, which are comparable with experimental ob¬ 
servations, are not explicitly related to fluctuations of 
the area enclosed by the two paths of the interferom¬ 
eter, because we used the same values of A l and A 
for all edge states composing the WP. 

As we can see from the model, in the mono-energetic 
plane-wave limit (a —)> oo) the energy-dependent ef¬ 
fects become less important, and we recover the ideal 
behavior with To = T\ = \ (notice that r(fco) = 
t(ko) ~ ^) and unitary visibility. The Gaussian 
modulation of the AB oscillations of the transmission 
coefficient disappears, since \F(k)\ 2 (27r)~ 1 S(k — 

ko): This is consistent with scattering-states models 
used in the literature 10,12,17 . 



Figure 6. (Color online) Transmission coefficient T cal¬ 
culated from the theoretical model for a = 20 nm, when 
scattering amplitudes t = r — ~^= are energy-independent 
(Eq. (25), bottom curve) and when scattering ampli¬ 
tudes are energy-dependent (Eq. (24), top curve). The 
energy-dependent curve has k e = 4.80 nm -1 , To = 0.757, 
Ti =s 0.245 and E = 40.4 nm, which are very close to the 
results of numerical simulations (see Table I). Note how 
the energy-dependent analytical model is able to repro¬ 
duce the numerical results of Fig. 3 (top curve) 


VI. CONCLUSIONS 

In conclusion, our time-dependent simulations al¬ 
low us to reproduce the interference pattern of elec¬ 
trons transmitted through Landau edge-states defin¬ 
ing a MZI. Specifically, we first used the split-step 
Fourier method to solve exactly the time-dependent 
Schrodinger equation, and then we calculated the 
single-particle dynamics and the transmission coeffi¬ 
cient of the interferometer, which shows Aharonov- 
Bohm oscillations. We exposed the effects of the fi¬ 
nite size of the electronic WP. In fact, the longitudi¬ 
nal extension (i.e. along the edge state) of the WP is 
connected to its energy uncertainty which, in turn, is 
determined by the injection process of the carriers: the 
larger the energy uncertainty, the more localized the 
WP. In Ref. 10, for example, the carrier injection in 
the device is achieved by an ohmic contact operating 
as the source terminal. However, the specific dynam¬ 
ics of the electrons moving from the Fermi sea of the 
metallic contact to the semiconductor edge state is un¬ 
known, in general. Indeed it is strongly affected by the 
atomistic details of the junction, by the temperature 
(i.e. the broadening of the Fermi level) and, in the 
case of Schottky contact, by the metal-semiconductor 
barrier. Furthermore, if the injection of the single 
carrier is realized via cyclic voltage pulses 25,53,54 , its 
WF can be partially tailored by the pulses timing and 
shape, and the correlation of several carriers, possi¬ 
bly levitons, can be estimated as reported in Ref. 55, 
thus assessing the transition from single-particle to 
multi-particle regimes.. With our approach, that in¬ 
cludes the energy selectivity of the QPCs, the real- 
space spreading of the carriers can be estimated a 
posteriori from the interference pattern 32,56 . 

In addition to simulating the exact dynamics, we de¬ 
veloped an analytical model that incorporates the 


















physics needed to reproduce the main features of the 
transport spectrum. Such a model must go beyond the 
standard delocalized scattering-state approach. In¬ 
deed, we derived Eqs. (24) and (25) by including the 
effects of the finite size of the electronic WP. Fur¬ 
thermore, the values of T 0 and T\ in Eq. (25) must 
depend on the different wave vectors fc composing the 
WP, in order to account for the energy-dependence of 
the scattering at the QPCs and, finally, to reproduce 
the offset value of the transmission, as a function of 
the spatial dispersion a. Also, the model leading to 
Eq. (24) can justify the reduction of the visibility of 
the interferometer with respect to the ideal case. 
Finally, we stress that the understanding of the dy¬ 
namics of single electrons in Landau ESs, possibly in¬ 
cluding 4-way splittings in QPCs and a non-trivial ge¬ 
ometry, is of utmost importance for several proposals 
of quantum computing architectures based on edge- 
channel transport. Usually, in these proposals, the 
quantum bits are encoded into two ESs of the same LL 
that are physically separated in space. However, some 
recent studies also focused on the possibility of cou¬ 
pling ESs belonging to different LLs 57-59 , between ESs 
belonging to two different 2DEGs through tunneling 60 
or between spin-resolved ESs of the same LL, which 
could also be used in the future to realize a different 
type of MZI 61,62 or scalable quantum gates. Also in 
the above systems, the effects of the spatial localiza¬ 
tion of the WP are pivotal for a realistic modeling of 
the devices. 
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it holds: 


1 


exp(<a(fc — fco)) + 1 


~ exp 


(a(k — fco) + e 


\2 1 


4e 


(A3) 


(here e is the Euler number). Therefore, since Eq. (24) 
for n=r 2 =r and t\ — t 2 = t can be written as 


T = 2tt 


J dk\F{k)f 


(A4) 


x ['r(fc) 4 + t(k ) 4 + 2r(fc) 2 t(fc) 2 cos(—fcA/ + A0)] , 


we deduce that, from each term in the square brack¬ 
ets, we can factor out exp(—2fy(& — fco) 2 ), which can 
be multiplied with |X(fc)| 2 oc exp(— 2cr 2 (k — fco) 2 ) to 
give an “effective” Gaussian weight with an increased 
standard deviation. We conclude that the energy- 
dependent value for E is given by: 


V 2 rv ^r 2 I 

E ~ a + —. 

4e 


(A5) 


With this approximated model we obtain, for a = 
20 nm, E = 39.7 nm, which is in excellent agreement 
with the value E = 40.98 nm obtained from the nu¬ 
merical fitting of the exact simulations (see Table I). 
A similar good agreement is obtained for a = 40 nm 
and a = 60 nm. 

Concerning the calculation of To and XI, we see from 
Eq. (A4) that the oscillating part of T and the Gaus¬ 
sian damping are given by the third term in the square 
brackets. Therefore, for large values of Al = 2a l, the 
rapid oscillations of the cosine cancel out, and only the 
first two terms survive, to give the saturation value of 
the transmission. Then 


The simulations were performed at the facilities of 
CINECA within the Iscra C projects “ESTENDI” 
and “TESSENA”. We would like to thank Arianna 
Bergonzini, Carlo Maria Bertoni, Carlo Jacoboni and 
Enrico Piccinini for useful discussions. 


Appendix A: Calculation of E, T 0 and XI 

In order to give an estimate for the energy-dependent 
value of E, we will use a perfectly symmetrical form 
for r(k) and t(fc), since the QPCs are half-reflecting 
for fc = fco and r(fc) 2 + t(k) 2 must be 1 for any fc; 
furthermore, we do not distinguish between QPC 1 
and 2: 

r(fc) = [exp(— a(k — fco)) + 1] _5 (Al) 

t(k) = [exp(+<a(fc - fc 0 )) + l] - ^ (A2) 

(where a = 1.13 • 10 _7 ra is obtained by fitting the nu¬ 
merical simulations). We now replace these functions 
with Gaussian profiles, since near the turning point 
fco and in the decaying tail of the Fermi distribution 


To = lim T (A6) 

a l —yoo 

= 27 t J dk |X(fc)| 2 [r(fc) 4 + £(fc) 4 ] . 

As we did for To, we can estimate XI by comparing 
Eq. (A4) with Eq. (25). As a first approximation, T\ is 
the maximum amplitude of the oscillating part, which 
is obtained when all the oscillating functions are max¬ 
imized and the Gaussian damping is zero, i.e. when 
the argument Acp of all the cosines in the integral van¬ 
ishes (and also Al = 0). Therefore, in this limit, we 
can subtract the background T 0 from T to get XI: 

XI = lim f - f 0 (A7) 

Ac /?—>-0 

= 2tt J dk\F(k )\ 2 [2 r(k) 2 t(k) 2 ] . 

For example, by using the above formulas, we obtain 
X 0 = 0.760 and XI = 0.240 for the WP with a = 
20 nm, that are in good agreement wit the values of 
Table I. It is straightforward to generalize these results 
to the cases when rq ^ r 2 and t\ ^ t 2 , to get the 
expressions of Eqs. (26) and (27). 
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SUPPLEMENTARY DATA 


In this Supplementary text, we show, as an example, the dynamics of a Lorentzian wave packet crossing a QPC 
of the kind described into Sec. II of the main text, and compare it with the corresponding result for a Gaussian 
wave packet, as used in the simulations of our work. 

We consider a ID Lorentzian wave function 




1 

(y - yo ) 



vt + isgn(kp)^ 


e ik F (y+vt ) 


where T is the full-width at middle height of the probability density |^f D (i/, t)| 2 , yo is the position of the 
maximum at the initial time t = 0. For a system with a linear dispersion E(k) = v • k, this corresponds to the 
wave function of a leviton 28 , but it describes also an electron emitted by a single-electron source quantum dot 
when subject to a suitable bias pulse 63 . The imaginary exponential kinetic factor exp (ikp(y + vt)) provides 
the group velocity (v g ) corresponding to the Fermi wave vector k F , which is given below. The factor sgn (k F ) 
is inserted to select only the wavevectors k associated with an energy E(k) > Ep , as required by the leviton 
energy profile. 

As in Eq. (5) of the main text, we Fourier transform in the k domain in order to obtain proper weights for the 
ESs k)e lky 


pL (k) = ^j dye- ik ytf D {y, 0 ) 

= _ ] j^^i e -i(k-k F )( 2 iyo+T Sg n(k F ))Q^ k _ k F ) S gn(k F )]sgn(k F ), 

where we used the Heaviside function ©[&]. Therefore, the expectation value for the group velocity is: 

(v 9 ) = h(k) = k F + sgn (k F )^, 


(SI) 


from which we can define the central wavevector ko = m ( v g ) /h of the wavepacket. With the approximation 
that all the involved edge states have the same transverse profile (which we used also in Sect. V to derive the 
analytical expression for the transmission coefficient), that is ip(x,k) ~ ip(x,k F ), the 2D wavepacket becomes 
(see also Eq. (4) of the main text) 


2p L (x, 


y) = f 


dkF L (k)<p hk (x)e l 


« <p{x, k F ) 



V™ (y ~Vo) +isgn(k F )^ 


(S2) 


thus showing that our linear superposition of ESs has the required Lorentzian shape. The above expression (S2) 
is taken as the initial Lorentzian wave packet and is propagated with the numerical approach described in the 
main text for the Gaussian case. Supplementary Figures 1 and 2 report the square modulus of the Lorentzian 
and Gaussian wave functions, respectively, at the initial time and at a following time step, after the action of 
the QPC. The wavepackets have been initialized to have the same group velocity (numerical parameters are 
given in the captions). Furthermore, Supplementary Figure 3 shows the effect of the QPC energy selectivity 
on the Lorentzian wave function. Although the time evolution in the latter case is more noisy than in the 
Gaussian case described in the main text, and the spatial spreading is clearly larger, the QPC is able to split 
the wave packet in two parts with essentially the same integrated probability of 0.5 but with different spectral 
components, as for the Gaussian case described in the main text. 
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Supplementary Figure 1. Square modulus of the Lorentzian wave function of Eq. (S2) at the initial time (a) and at 
t = 6ps, when the WP has been split in a reflected (b) and a transmitted (c) component by the QPC. Here Y = 40 nm 
and the wavevector ko associated with the group velocity is ~ 1.2 nm -1 , while the other parameters are the same as in 
Fig. 1 of the main text. 
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Supplementary Figure 2. Square modulus of the Gaussian wave function of Eq. 4 of the main text at the initial time (a) 
and at t — 5ps, when the WP has been split in a reflected (b) and a transmitted (c) component by the QPC. We used 
the same other parameters as in Fig. 1 of the main text (i.e., a = 40 nm and ko ~ 1.2 nm -1 ). 
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Supplementary Figure 3. Reciprocal-space weights, as in Eq. (SI), of the Lorentzian wave function at the initial time 
(yellow) and for the transmitted (blue) and reflected (red) components after the split at the QPC. Their real-space 
evolution is shown in Supplementary Figure 1. 































